Homogeneous nucleation of a non-critical phase near a continuous phase transition 
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' Homogeneous nucleation of a new phase near a second, continuous, transition, is considered. The 

continuous transition is in the metastable region associated with the first-order phase transition, 
one of whose coexisting phases is nucleating. Mean-field calculations show that as the continuous 
transition is approached, the size of the nucleus varies as the response function of the order parameter 
of the continuous transition. This response function diverges at the continuous transition, as does 
the temperature derivative of the free energy barrier to nucleation. This rapid drop of the barrier 
as the continuous transition is approached means that the continuous transition acts to reduce the 
barrier to nucleation at the first-order transition. This may be useful in the crystallisation of globular 
i-C • proteins. 

o 

I. INTRODUCTION 

> 

■*_> | When a phase transition is first order the formation of a new phase is an activated process [jlj . A nucleus of the new 
phase must form, and overcome a free-energy barrier before it can grow into a new phase. The rate at which nuclei 
overcome a barrier of height AF* scales as exp(— AF*/kT) and so this rate is very sensitive to the barrier height 
. For definiteness consider a first-order phase transition at a temperature T a in which the low temperature phase 
is the ordered phase, and the high temperature phase is the disordered phase. If we cool the disordered phase below 
T a but the barrier to nucleation of the ordered phase is very high then the rate at which nuclei of the ordered phase 
£h ' form will be effectively zero and the ordered phase will not form even though its free energy is lower than that of the 
disordered phase. The disordered phase will persist; it is said to be metastable. Here we calculate AF* for nucleation 
i^i - near a second, continuous, transition, which we call transition j3. Continuous transitions are critical points and so 
exhibit universal and beautiful behaviour, the thermodynamic and correlation functions contain power law terms with 
t-H , exponents which depend only on dimensionality and the symmetry of the order parameter ||. A priori, we might 
expect that AF* might also contain a power law term with an exponent which depends only on dimensionality and 
the symmetry of the order parameter. This would allow us to make predictions about how AF* varied near a critical 
point which would apply to a whole class of systems. Below, we present the results of calculations within mean-field 
theory, for an Ising-like continuous transition in three dimensions. We determine the singular power law part of the 
free energy barrier AF*: just above Tp, it varies as where £ is the correlation length associated with the order 
parameter of transition j3. This singular part means that the derivative of AF* with respect to temperature diverges 
as the critical point is approached: the barrier to nucleation drops rapidly just above the transition. This agrees with 
the pioneering simulations of ten Wolde and Frenkel j^] who found an anomalously low AF* near the critical point 
5 , of a metastable fluid-fluid transition. 

As far as the author is aware, nucleation near a critical point has only been considered for a fluid-fluid critical point 
within a strongly first-order fluid-crystal transition. This work was inspired by the observation that globular proteins 
often crystallise at temperatures close to where we might expect a metastable fluid-fluid critical point |6[ [?J , and that 
the phase diagrams of some globular proteins do possess a metastable fluid- fluid transition ||,||]. But as the effect we 
. will examine is due to the decreasing free energy cost of fluctuations near a critical point it is universal, i.e., applies 
to any other system in the same universality class, that of the three-dimensional Ising model. Indeed, it is easy to 
show that it also applies to systems in the universality classes of the Ising model in other dimensionalities . See 

■ Refs. [pr|-Jl5[ and references therein for recent work. 

The next section briefly sets out the standard Landau theory for a continuous transition with a scalar order 
parameter. Section III then calculates the free energy of a nucleus within a simple mean-field theory of the square- 
gradient type. Derivatives with respect to temperature and external field are also found. The final section is a 
conclusion. See Refs. |||4|] for an introduction to systems near critical points and Ref. M for an introduction to 
homogeneous nucleation. 
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II. BULK BEHAVIOUR 



We have a system, which at equilibrium has one phase transition: a strongly first order transition, transition a, 
which is at a temperature T a . For definiteness we let the high-temperature phase be the disordered phase and the 
low-temperature phase be the ordered phase. If we consider a very pure sample JO] then we can supercool down to 
temperatures below T a to obtain metastable states 0, i.e., the disordered phase is stable for long (with respect to the 
relaxation time of the system) periods of time over some temperature range just below T a . It is stable because the 
formation of the ordered phase is an activated process, the ordered phase must nucleate, overcoming some free energy 
barrier AF* which will be a strong function of temperature and which diverges as T — > T~. Here, we are interested 
in how AF* behaves if as we cool down in the metastable disordered phase, we approach a continuous transition, 
transition /3, at some temperature Tp <T a . 

We will assume that the nucleation barrier to transition a is so large that it is possible to start from some temperature 
T > T a and slowly cool down past T a either to and below Tp, or at least to a temperature only a little above Tp, 
without transition a occurring. If it is possible to cool slowly down to Tp, then transition (3 is said to be metastable [0: 
it is observable. If it is not possible to cool slowly down to Tp without transition a occurring then clearly transition 
f3 is not observable; it is unstable not metastable with respect to transition a ||. We will be studying nucleation of 
the low temperature phase of transition a as Tp is approached from above and so we will not only be determining 
the effect of transition (3 on a but also looking at whether or not (3 is observable. Roughly speaking, if the proximity 
of a transition (3 acts to strongly enhance the nucleation rate of transition a, then this nucleation rate may become 
large thus rendering transition (3 unobservable. We should also mention that we are using temperature as our variable 
simply for definiteness, we could replace it by another field variable, e.g., pressure, without changing our conclusions. 

So, starting from high temperature and cooling down below T a we can approach Tp. The order parameter of 
transition (3 is denoted by to and it may or may not be related to that of transition a. The external field which 
couples to m is h. The theory here will be mean-field in nature but rather general. We only have to assume that 
the nucleus of the ordered phase of transition a has a core which has properties close to that of the bulk ordered 
phase (c.f., the assumptions which underly classical nucleation theory 0) and that this core couples to the order 
parameter of transition (3. By coupling to m we mean that if there is a nucleus at the origin, then the local value 
of m, m r (r) ^ m, where r is the distance from the centre of the nucleus. Both assumptions are very reasonable: 
for a strongly first-order phase transition it is difficult to imagine a situation where the nucleus does not have a core 
with near bulk properties, and the core of the nucleus must perturb its surroundings and so, in the absence of special 
symmetries, will locally perturb the order parameter of transition (3. Figure [l] is a schematic of the nucleus. 

Near Tp we use a Landau theory for the transition (3. The Landau theory of a continuous transition is simple, it 
is a textbook problem, see for example Chaikin and Lubensky's j| or Kadanoff's j|. The bulk free energy per unit 
volume /(to) as a function of the order parameter m is 

/(m) = ^a(T - Tp)m 2 + m 4 - hm, (1) 

The transition is at Tp at h = 0. We will only examine behaviour at h = but we retain h in order to look at the 
response of the system to an external field which couples to to. Below, when we study the nucleus near transition (3, 
we will find that the outermost part of the density profile of the nucleus is controlled by the response function of to, 
X, defined by 

_ x = (3h^ = 

which is, using Eq. (Il 



' dm J \ dm 2 ) ' 



X- 1 = a(T -Tp) T > Tp. (3) 



III. THE NUCLEUS 



We split the nucleus into two parts: the core and the fringe. The core is that part of the nucleus less than r c from 
its centre and the fringe is that part farther away than r c . The fringe of the nucleus is assumed to be spherically 
symmetric. The core of the nucleus contains at its centre a volume which is close to the bulk ordered phase of 
transition a. The fringe is essentially the region which surrounds this core and is perturbed by the core. Its radius 
is therefore the correlation length £ of to and so diverges as Tp is approached. As we are concentrating on universal 
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aspects of the nucleus and of AF* we will replace the core by a boundary condition on m r in the fringe. We set 
m r (r = r c ) = m c which is taken to be independent of temperature and of h. Also, m r (r — > oo) = m, which is just 
the obvious boundary condition that to must tend towards its bulk value far from the nucleus. Note that as we are 
above the transition and are working at zero field, to — but we will retain an explicit m dependence in order to be 
able to take derivatives of AF*. In the fringe and near Tp, we need only consider the order parameter for transition 
(3 and the variations in m r will not be large. Therefore, we employ a standard Landau-Ginsburg or square-gradient 
functional for the excess free energy AF of an inhomogeneity in an otherwise homogeneous phase, P, [lq [l8[ 



AF = AF r 



Af(m r ) + k (S7m r y 



dr, 



where 



Af(m r ) = f(m r ) - /(to) - h(m r - to), 



(4) 



(5) 



is the bulk free energy change per unit volume to go from to to m r . The excess free energy AF is the free energy with 
a nucleus minus that without a nucleus; AF C is the contribution of the core. The second term within the brackets of 
Eq. (Q) is a gradient term: the free energy cost due to variations in space of m r . It is proportional to the gradient 
squared which is the lowest order term in a gradient expansion and so is only adequate when to is slowly varying. The 
coefficient, k, of this term is taken to be a constant. The total excess of m due to a nucleus is equal to the integrated 
value of m r — to. Defining Am(r) — m r (r) — to, then the integral over all space of this function gives the total excess 
of the order parameter due to the nucleus, 

Ato = J Am(r)dr. (6) 

The free energy barrier AF* is the value of AF for the nucleus when it is at its maximum, at the top of the 
barrier. The nucleus at the top of the barrier is called the critical nucleus p9|| . For the critical nucleus we may set 
the functional derivative of AF with respect to the profile m r (r) to zero, 



fdAf(r 



2K\7 2 m r 



r > r c 



(J) 



Once we have solved Eq. (fy we can insert the solution into Eq. (Q) to obtain the excess free energy of the critical 
nucleus, AF*. 

The fringe is the outermost part of the nucleus, where m r is near the bulk value to. So we can use a Taylor 
expansion about m r = m for A/, 



A/K) = ^ (Ato) 2 



dAf(m r ) 
dm 



= x 1 Am 



because both A/ and its first derivative are zero for m r = to, and the second derivative is x 1 [Eq. 
Eq. (|) into Eq. © we have 

X -1 Am(r) - 2kV 2 Ato(t-) = 0, 
which has a solution of the Ornstein-Zernike form, 

Am(r) = (to c - to) ^ exp [(r c - r)/£] , 

with £ the correlation length for to, given by 

t = (2*X) 1/2 

£= (2K/a) 1/2 (T-T p y 1/2 T>T p , 



(8) 
(9) 

Substituting 
(10) 

(11) 



(12) 
(13) 



where Eq. (12) defines £ and we used Eq. (||) to obtain expressions for £ near Tp, Eq. @. To obtain Eq. (|ll|) the 
boundary conditions m r {r — ► oo) = to and m(r c ) = m c were employed. It is not necessary to specify r c or m c beyond 
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saying that they should be such that m c — m is small and so Am(r) will, as required for Eq. jic|), be small for r > r c . 
From Eq. (|ll| ) we see that the width of the fringe is, as we expected, of the order of the correlation length £ for m. 

Having obtained the density profile, Eq. (O), we can substitute this into Eq. (0), using Eq. (||) for A/, and obtain 
an expression for the free energy barrier to nuclcation. We have 



poo 

AF* = AF C + Anr 2 c (m c - m) 2 J dr 

poo 

= AF C + Airnr 2 (m c - m) 2 dr 



2 X 
2 

e 

AF C + 4nKr c (m c - m) 2 [1 + r c /£] , 



1 1 

r i 



2 



1 



exp [2(r c - r)/i] 
exp[2(r c -r)/£] 



(14) 



where in obtaining the second line from the first we substituted £ for x using Eq. (|12|). Finally, we can set m = to 
obtain the free energy barrier to nucleation of the ordered phase of transition a near transition f3, 



AF* = AF C + 4irKr c m 2 c [1 + r c /£] . 
As we approach transition (3, T — > Tp, AF* approaches the finite limit 

AF*(T ^T ) = AF C + 4TTnr c m 2 c . 

The free energy AF* can be written as 

AF* = AF*(T = T (j )+A 



(15) 



(16) 



(17) 



where A = 4:irKr c m 2 , a constant. The singular part of A_F* has the form: the ratio r c /£ raised to the power 1. 



A. Derivatives of AF* 

We can take the temperature derivative of AF* . As m does not vary with T above Tp we may use Eq. ([l5]) , and 
obtain 

OAF* 8AF 2 2 flr x n «s 

which near Tp becomes 

A (AF* - AF C ) = {2na) 1/2 irr 2 m 2 c (T - T p y 1/2 T > Tp, (19) 

where we used Eq. ( |l3|) for £. Just above the transition j3 the derivative of the barrier diverges to +oo; the barrier 
drops very rapidly with decreasing temperature just above Tp. 

We can also take the derivative of AF* with respect to the field h conjugate to m. Using Eq. (|l4]), and taking note 
of the definition of x-, Eq- ([§), 

dAF* dAF r 

—r^- = - STT K r c m cX [1 + r c /£] , (20) 

where after taking the derivative we set m = 0. Note that d^/dh — 0. As transition is approached the rate of 
change of AF* with respect to the field conjugate to the order parameter diverges as the response function x- Also, 
if we substitute our solution for m(r > r c ), Eq. (|ll]), with m — 0, into Eq. (|^), we obtain the size of the nucleus 

Am* = Am* + %-Knr c m c x [1 + r c /£] (21) 

with 

Am c = / Am(r)dr, (22) 
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the contribution of the core, and we used Eq. ( |l2| ) to substitute x f° r £ 2 - Comparing Eqs. (20) and (pl|), we see that 

JL (AF* - AF C ) = - (Am* - Am c ) . (23) 

For the fringe, the derivative of the free energy with respect to h is equal to minus the excess m. This result is 
essentially what is called the nucleation theorem |2(] [2^] in studies of nucleation in fluids. It states that the larger 
the nucleus, the larger Am* is, the more rapidly the nucleation barrier varies with h. In fluids m is a number density 
difference and h is the chemical potential. 

Returning to Eq. (|2l]) for Am* we see that although the core can only contribute a finite amount to Am* as its 
volume is finite, the contribution of the fringe diverges as transition (3 is approached. The size of the nucleus diverges 
as x as the continuous transition is approached. This result was first derived by the author in Ref. fll5| ]. See also the 
earlier work of Talanquer and Oxtoby jllj who first suggested that the size of the nucleus diverges as a critical point 
is approached. In Refs. the critical point is that of a fluid-fluid or vapour-liquid-like transition. 



IV. CONCLUSION 



We have considered the effect of a continuous transition, transition /3, on the homogeneous nucleation of a new 
phase at a first-order transition, transition a. We found that the temperature derivative of the free energy barrier to 
nucleation, AF*, diverged as (T — Tp)^ 1 / 2 within our mean-field theory, and that the size of the critical nucleus, the 
nucleus at the top of the free energy barrier to nucleation, diverged as the response function \ ~ (T — Tp)~ l , The 
presence of a critical point makes the nucleus very large, its diameter is the correlation length £, and causes the free 
energy barrier to nucleation to decrease rapidly with decreasing temperature. It reduces the barrier and so facilitates 
nucleation. This is just what was first demonstrated by ten Wolde and Frenkel || for nucleation of a crystalline phase 
near the critical point of a fluid-fluid transition. It is a rather general phenomenon and applies to any continuous 
transition with a scalar order parameter, i.e., any Ising-like transition. Whether or not the same effect appears near 
a continuous transition in a system which is anisotropic or in which the order parameter is not a simple scalar, is an 
interesting open question. 

It is a pleasure to acknowledge discussions with R. Evans and A. Parry. Work supported by EPSRC (GR/N36981). 



Here we consider only homogeneous nucleation, nucleation in the bulk, far from any interface. Heterogeneous nucleation, 
nucleation at an interface, either with a wall or an impurity, is in fact more common, and is not always an activated 
process. 

P. G. Debenedetti, Metastable Liquids (Princeton University Press, Princeton, 1996). 
L. P. Kadanoff, Statistical Physics (World Scientific, Singapore, 2000). 

P. M. Chaikin and T. C. Lubensky, Principles of Condensed Matter Physics (Cambridge University Press, Cambridge, 
1995). 

P. R. ten Wolde and D. Frenkel, Science 277, 1975 (1997). 

A. George and W. Wilson, Acta. Crystallogr. D 50, 361 (1994). 

D. Rosenbaum, P. C. Zamora and C. F. Zukoski, Phys. Rev. Lett. 76, 150 (1996). 

M. L. Broide, C. R. Berland, J. Pande, O. O. Ogun and G. B. Benedek, Proc. Nat. Acad. Sci. 88, 5660 (1991). 
M. Muschol and F. Rosenberger, J. Chem. Phys. 107, 1953 (1997). 
R. P. Sear, unpublished work. 

V. Talanquer and D. W. Oxtoby, J. Chem. Phys. 109, 223 (1998). 

C. Haas and J. Drenth, J. Phys. Chem. B 104, 368 (2000). 

N. M. Dixit and C. F. Zukoski, J. Coll. Int. Sci. 228, 359 (2000). 

D. Pini, G. Jialin, A. Parola and L. Reatto, Chem. Phy s. Lett. 327, 209 (2000). 



R. P. Sear, J. Chem. Phys, in press; cond-mat/0010154 



J. W. Cahn and J. E. Hilliard, J. Chem. Phys. 31, 688 (1959). 
J. W. Cahn and J. E. Hilliard, J. Chem. Phys. 28, 258 (1958). 
B. Widom, Chem. Soc. Rev. 14, 121 (1985). 

Note that conventionally 'critical' is used to denote both a continuous transition and the crystallite at the top of the barrier 
although there is no connection between the two uses of the word. Although this terminology is a little unfortunate we 
will use it here and so we will be studying a critical nucleus near a critical point. 



5 



[20] D. Kashchiev, J. Chem. Phys. 76, 5098 (1982). 

[21] Y. Viisanen, R. Strey and H. Reiss, J. Chem. Phys. 99, 4680 (1993). 

[22] R. K. Bowles, R. McGraw, P. Schaaf, B. Senger, J.-C. Voegel and H. Reiss, J. Chem. Phys. 113, 4524 (2000). 



FIG. 1. Schematic of a nucleus of the ordered phase of transition a near transition /3. The core of the ordered phase of 
transition a is solid black, and the perturbation this causes in the surroundings is the shaded circle of radius the correlation 
length £. The sphere with radius r c which divides the nucleus into a core and surroundings is denoted by a dashed circle. 
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